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: 

■ The topic of this conference is "The Chaotic Universe". One of the main achievements of last 
^ ' century has been to relate chaos in fluids to their thermodynamics. It is our purpose to make 
Q , connection between chaos in gravitation and standard thermodynamics. Though there have been 

•"7" ■ many previous steps and attempts, so far no convincing conclusion has been reached. 

' After explaining how the approach works for glasses, we shall discuss the thermodynamics of 

0^ , two speciflc systems: black holes and globular star clusters. In both cases we point out that the 

■ dynamics satisfies the first and second law of thermodynamics, though negative specific heats occur. 

o ■ 

CD ■ I. INTRODUCTION 

^ _ Thermodynamics is the old science that describes the flow of energy in systems with many atoms. Works started 
. by Carnot, Kelvin and Clausius showed that these laws are very general. This universality led to the formulation of 
C/3 ' the first and second law of thermodynamics, that apply to a vast amount of systems, such as gases and crystals. 

During half a century there was still a problem with the application to glasses. In this field there were classical 
paradoxes related to the so-called Ehrenfest relations and Prigogine-Defay ratio. The solution of this problem is 
discussed below. Due to its inherent non-equilibrium nature a glass is far from equilibrium. To describe it in 

■ a thermodynamic treatment one has to take into account at least one additional system parameter, the effective 
temperature, and its conjugate variable, the configurational entropy |]l| 

Q After having realized how thermodynamics should be formulated for glasses, we have investigated the situation for 

^ black holes . For this problem various aspects of the dynamics are known, and it was generally expected that the 
laws for black hole dynamics would coincide with the laws of black hole thermodynamics. We nevertheless felt that 
I ■ the proper connection between black hole dynamics and standard thermodynamics had not been made, and this will 
^ , be clarified below. 

^-H ' The question whether thermodynamics also applies to other self-gravitating systems, such as star clusters, is still 
' open today. A major fraction of researchers believes that the final answer will be negative. Their reservations are 
partly based on non-extensivity of thermodynamic variables due to the long range gravitational potential, and the 
divergence of the partition sum due to the strong short-distance attraction, see Q for a review. We shall discuss, 

■ however, that the general formulation of non-equilibrium thermodynamics applies on timescales where evaporation 
Qs^ occurs, and partition sums cannot be used. 

II. THERMODYNAMICS FAR FROM EQUILIBRIUM 

' ^ , In the gravitation literature the first law of thermodynamics is usually presented as dU = TdS — pdV, and the 
Q ■ second law as dS' > 0. These non-general formulations apply only to equilibrium and to a closed system, respectively. 
^ ' We should stress that they do not apply to many situations far from equilibrium. This has often led to the belief that 
thermodynamics does not work. However, to investigate that properly, we must investigate whether the most general 
formulation does apply. 

. The general formulation of the first law of thermodynamics merely says that energy is conserved: the change of the 
$H ' systems energy equals the heat added to the system plus the work done on the system, i.e. 

a ! 

dU =dQ+dW (1) 

The second law says that disorder can only increase. Heat can only go from high temperatures to low temperatures, 
which implies 

dQ < TdS (2) 
where S is the total entropy. The equality sign holds if and only if there is equilibrium. 
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III. GLASSES 



Thermodynamics for systems far from equilibrium has long been a field of confusion. A typical application is 
window glass. Such a system is far from equilibrium: a cubic micron of glass is neither a crystal nor an ordinary 
under-cooled liquid. It is an under-cooled liquid that, in the glass formation process, has fallen out of its meta-stable 
liquid equilibrium. There is thus a separation of timescales between the fast processes, often called /3-processes, and 
the basicaly quenched a- or configurational processes. 

Until our recent works on this field, the general consensus reached after more than half a century of research was: 
Thermodynamics does not work for glasses, because there is no equilibrium This conclusion was mainly based on 
the failure to understand the Ehrcnfest relations and the related Prigogine-Defay ratio. It should be kept in mind 
that, so far, the approaches leaned very much on equilibrium ideas, q], |@] [||. We shall stress that such approaches 
are not applicable, due to the inherent non-equilibrium character of the glassy state. Even the quoted conclusion 
itself is actually confusing, since thermodynamics should also hold outside equilibrium. 

Thermodynamics is the most robust field of physics. Its failure to describe the glassy state is quite unsatisfac- 
tory, since up to 25 decades in time can be involved. Naively we expect that each decade has its own dynamics, 
basically independent of the other ones. We have found support for this point in models that can be solved exactly. 
Thermodynamics then means a description of system properties under smooth enough non-equilibrium conditions. 



A state that slowly relaxes to equilibrium is characterized by the time elapsed so far, sometimes called "age" or 
"waiting time" . For glassy systems this is of special relevance. We shall restrict to systems with one diverging time 
scale. They are systems with first-order-type transitions, with (smoothed out) discontinuous order parameter, though 
usually there is no latent heat. 

The picture to be investigated in this work starts by describing a non-equilibrium state characterized by three 
parameters, namely T,p and the age t. By exactly solving the dynamics of certain model systems [ pT| , 
it was found that in the thermodynamics the non-equilibrium nature shows up only through an effective temperature 
Te{t); also numerical data on the slow evolution of a binary Lennard- Jones glass were recently interpreted in terms of an 
effective temperature . For a set of smoothly related cooling experiments Ti (t) at pressures pi , one may express the 
effective temperature as a continuous function: Te,i{t) Te{T,p). This sets a surface in (T, Te,p) space, that becomes 
multi- valued if one first cools, and then heats. For covering the whole space one needs to do many experiments, e.g., 
at different pressures and different cooling rates. The results should agree with findings from heating experiments and 
aging experiments. Thermodynamics amounts to giving differential relations between observables at nearby points in 
this space. In principle also an effective pressure could be needed. If there are many long time scales, also several 
effective temperatures would be needed. We shall restrict ourselves to the simplest case of one effective temperature. 

Of special importance is the thermodynamics of a thermal body at temperature T2 in a heat bath at temperature 
Ti . This could apply to mundane situations such as a cup of coffee, or an ice-cream, in a room. There are also two 
entropies, Si and 82- Notice that there are also two time-scales: the time-scale for heat to leave the cup is much 
larger than the time-scale for equilibrating that heat in the room. It is this separation of time-scales that makes 
spontaneous difference in temperatures possible. The change in heat of such a system obeys dQ < TidSi + T2dS2- 

A similar two-temperature approach proves to be relevant for glassy systems. The known exact results on the 
thermodynamics of systems without currents can be summarized by the very same change in heat H |lt| 



where Sep is the entropy of the fast or equilibrium processes (/3-processes) and T the configurational entropy of the 
slow or configurational processes (a-processes). This object is also known as information entropy or complexity. 
Notice that it deviates from the standard definition j^] that the configurational entropy Sc is the entropy of the glass 
minus the one of the vibrational modes of the crystal. The latter definition also incorporates fast non-vibrational 
processes. 

It is both surprising and satisfactory that a glass can be described by the same general law. If, in certain systems, 
also an effective pressure or field would be needed, then dQ is expected to keep the same form, but dW would change 
from its standard value —pdV for liquids. Such an extension could be needed to describe a larger class of systems. 

For a glass forming liquid the first law dU — dQ + dW becomes 



A. Two-temperature thermodynamics 



dQ = Td5op + TedJ 



(3) 



dU = Td5cp + Tedl - pdV 



(4) 
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It is appropriate to define the generalized free enthalpy 

G = U - TS^p - TJ + pV 
This is not the standard form, since ^ T. It satisfies 

dG = -S'epdT - IdTe + Vdp 

The total entropy is 



The second law requires dQ < TdS, leading to 



S — Sep -\- X 



(Te - T)dX < 0, 



which merely says that heat goes from high to low temperatures. 

Since — Te{T,p), and both entropies arc functions of T, and p, the expression yields a specific heat 



dQ 
dT 



p ^ dT 



T,p dT 



dl 

T^.P^dfe 



T,p dT 



dT, 

In the glass transition region all factors, except drT^, are basically constant. This leads to 

dT 



) 



Cp — Ci + C2 



dT 



(5) 
(6) 
(7) 
(8) 

(9) 
(10) 



Precisely this form has been assumed half a century ago by Tool fl^ as starting point for the study of caloric behavior 
in the glass formation region, and has often been used for the explanation of experiments Q p4[ . It is a direct 
consequence of eq. (0) . 



B. Modified Maxwell relation, Ehrenfest relation and Prigogine-Defay ratio 



For a smooth sequence of cooling procedures of a glassy liquid, eq. (jj) implies a modified Maxwell relation between 
macroscopic observables such as U{t,p) U{T,p) = U{T,Te{T,p),p) and V. This solely occurs since Tg is a non- 
trivial function of T and p for the smooth set of experiments under consideration. Most Maxwell relations involve the 
entropy, which is difficult to observe. The modified Maxwell relation between the observables U and V reads 



dU 
dp 



dV 

T op 



T dT 



dT 



dT, 
p dp 



T op 



dT, 
T dT 



dp 



(11) 



In equilibrium T, — T, so the right hand side vanishes, and the standard form is recovered. 

In the glass transition region a glass forming liquid exhibits smeared jumps in the specific heat Cp = d{U+pV) /9T|p, 
the expansivity a = d\nV/ T\p and the compressibility k ~ ~d\nV/dp\T- If one forgets about the smearing, one 
may consider them as true discontinuities, yielding an analogy with continuous phase transitions of the classical type. 

Following Ehrenfest one may take the derivative of /Sy{T,pg{T)) = 0. The result is 



dT 



(12) 



The conclusion drawn from half a century of research on glass forming liquids is that this relation is never satisfied 
[H [|6| j^. This has very much hindered progress on a thermodynamical approach. However, from a theoretical 
viewpoint it is hard to imagine that something could go wrong when just taking a derivative. We have pointed out 
that this relation is indeed satisfied automatically but it is important say what is meant by k in the glassy state. 

Previous claims about the violation of the first Ehrenfest relation can be traced back to the equilibrium thermody- 
namics idea that there is one, ideal k, to be inserted in (p^). Indeed, investigators always considered cooling curves 
V{T,pi) at a set of pressures pi to determine Aa and dpg/dT. However, Ak was always determined in another way, 
often from measurements of the speed of sound. In equilibrium such alternative determinations would yield the same 
outcome. In glasses this is not the case: the speed of sound is a short-time process, which only measures short-time 
relaxation. Therefore alternative procedures should be avoided, and only the cooling curves V{T,pi) should be used. 
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They constitute a liquid surface Miquid p) and a glass surface Vg\ass{T,p) in {T,p, V) space. These surfaces intersect, 
and the first Ehrenfest relation is no more than a mathematical identity about the intersection line of these surfaces. 
It is therefore a tautology as was also stressed by McKenna |17|. 

The second Ehrenfest relation follows from differentiating AU{T,pg{T)) — 0. The obtained relation will also be 
satisfied automatically. However, one then eliminates dU /dp by means of the Maxwell relation (|ll|). We obtain 



TgV dT y V dT 



dX 
dT 



(13) 



The last term, involving the "total" derivative dX/dT ~ dl/dT + {dX/dp){dpg/dT) of the configurational entropy 
along the glass transition line, is new. Its prefactor only vanishes at equilibrium, in which case the standard Ehrenfest 
relation is recovered. The equality Te{T,pg{T)) = T implies the useful identity 



dTe 
dT 



dT 



+ 



dT, 
dp 



dpg_ 

T dT 



Historically one has defined Prigogine-Defay ratio 



^ ACpAK 



TV{Aa)^ 



(14) 



(15) 



This looks like an equilibrium quantity, and for equilibrium transitions it should be equal to unity. Assuming that at 
the glass transition a number of unspecified parameters undergo a phase transition, Davies and Jones derived that 
n > 1 while DiMarzio showed that in that case the correct value is H = 1 [|l^. In glasses typical experimental 
values are reported in the range 2 < H < 5. It was therefore generally expected that H > 1 is a strict inequality 
arising from the requirement of mechanical stability. 

However, since the first Ehrenfest relation is satisfied, it holds that H 



H 



ACr, 



TVAa{dpg/dT) 



1 



1 



VAc 



1 



dTe 
dT 



dl 
dp 



(16) 



Depending on the smooth set of experiments to be performed, dpg/dT can be small or large: H depends on the 
set of experiments. As a result, it can also be below unity. Rehage-Oels found H — 1.09 « 1 at p = 1 kbar 
using a short-time value for k. Reanalyzing their data we find from (p^, where the correct n has been inserted, a 
value H = 0.77, which indeed is below unity. The commonly accepted inequality H > 1 is based on the equilibrium 
assumption of a unique k. Our theoretical arguments and the Rehage-Oels data show that such assumptions are 
incorrect. 

Further steps involve fluctuation formula and the fluctuation-dissipation theorem. They are modified outside 
equilibrium, and effective temperatures also occur there. In sirnple model systems all these effective temperatures 
coincide, and this is expected to hold for a subclass of systems M |ll[]. 



IV. BLACK HOLES 



Black holes were anticipated by Laplace in 1798. He considered the gravitational escape problem from a mass 
M. Equating kinetic mv'^/2 to the potential energy GMm/r he observed, for light, the critical escape radius Rs = 
2GM/c?. Exactly this radius shows up in the Schwarzschild solution of the Einstein equations. The metric for a 
neutral, non-rotating black hole reads 



ds' 



1 



(17) 



1 - Rs/r 

0, implying a radial speed dr/dt = (1 — Rs/r) c, which vanishes at 



For spherical light waves one has ds = d6' = d 
the 'horizon' r = Rg. 

The connection between black holes and thermodynamics goes back to Bekenstein. He introduced in 1973 the 
notion of black hole entropy, proportional to its area, in dimensionless units |po| . Since Hawking later fixed the 
numerical prefactor, the result is now called the Bekenstein-Hawking entropy 



Sbh — 



fiG 



(18) 
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The presence of Ti calls for a quantum mechanical interpretation, and the quantum evaporation of black holes was 
demonstrated by Hawking This underlined the physical relevance of Bekenstein's approach. The black hole 

radiates as a black body at Hawking temperature 

TiGk hc'^ 

" ^ 2TTc^kB " SirGMkB ' 

where the second equality holds for a non-rotating, neutral black hole. All possible particles are emitted at this 
temperature; for large black holes, however, Th is so small, that in practice only massless particles (photons, gravitons, 
and possibly neutrino's) are emitted. 

A black hole has no hair, i.e. it can be characterized by a few parameters, namely its mass M, charge q and angular 
momentum J. This is reminiscent of fluids, that can be characterized by temperature and pressure. From the mass 
formula for charged, rotating black holes |2|] , it is known for long that the energy U = Mc^ satisfies 

dU ^ ^dA + Q-dJ + (l)dq (20) 
ovr 

where k is the surface tension, A = 47ri?| the area, the horizon's angular velocity, and cj) the electrostatic potential 
at the horizon. This law holds when adding matter to one given black hole, but also when comparing two different 
black holes. These two very different applications suggest a universal validity, and a thermodynamic description. 

These two fundamental results led Bardeen, Carter and Hawking |2^ to formulate "the four laws of black hole 
dynamics" . The zeroth law states that the surface tension k is constant at the black hole surface, just like the 
temperature is the same everywhere in an equilibrium system. The first law is given in eq. (|2l|). Since the last two 
terms corresponds to work terms, one may write this relation in the suggestive form 

dU ^TndSBH +dW (21) 

This formulation is sometimes referred to as the first law of black holes thermodynamics, but so far it is only a 
rewriting of (|20|). Bekenstein had already discussed the generalized second law of black hole dynamics: the total 
entropy S = Sbh + ^rm where 5',„ is the entropy of the matter outside the black hole, cannot decrease: 

d5 > (22) 

The third law concerns "extremal" black holes, the ones that are maximally rotating and/or are maximally charged, 
and have k = 0. The third law of black hole dynamics says that black holes with Th = cannot be reached by a 
finite number of steps |Q . 

As indicated by their name, the "four laws of black hole dynamics" look similar to the laws of thermodynamics, 
though a precise connection was not established. Only for a black hole put in a not-too-large box there is equilibrium, 
and the standard laws of equilibrium thermodynamics apply ||2^ Q . During last 25 years the outstanding question 
has been to relate these laws to the standard laws of thermodynamics. 

From the view point of a condensed matter physicist, the literature on black hole thermodynamics is somewhat 
confusing. First of all, one should define the system for which a thermodynamic description is to be given. A 
natural choice is to consider as system the black hole and a "Gedanken" sphere around it of, say, a hundred times 
the Schwarzschild radius. One could also consider the whole universe as an isolated container. Our next objection 
concerns the non-general formulation of the first and second law. This has already be discussed in a more general 
context above. 

Having defined the system, one should determine its entropy. For the Gedanken sphere with the black hole in it, eq. 
( ^ ) applies. The entropy occurring in eq. ( p2| ) is 5 = Sbh + S^^. The latter is the entropy of the cosmic background 
matter outside the black hole but inside the Gedanken sphere, and expected to scale with the sphere's volume. The 
radiation generated by the hole will quickly leave the system and go to the heat bath around it; this is described by 
adQ < 0. 

If, on the other hand, the whole universe is considered as system, then dQ = 0. If no work is done, this implies 
that dU ~ 0, saying that energy radiated from the hole is still inside the system. In that case eq. ( [2l| ) does not 
describe the change of the system's energy, it only says something about the black hole. The total entropy is now 
S = Sbh + Sm, and the second law of thermodynamics indeed says that dS > 0. 

We conclude that eq. (|2^) and ( p^ are both valid, though they should not be applied simultaneously. They refer 
to different situations, that is to say, to different time scales. When only the black hole and its Gedanken sphere 
are considered, this describes the radiation emitted in per unit time. When considering the change in entropy of 
the whole universe, one tacitly assumes time scales so large that all information on the emitted radiation has been 
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lost. Notice that the negative curvature of the universe, leading to exponential divergence of nearby trajectories, can 
establish this loss of information even though the emitted photons are hardly scattered. 

A final, severe, objection against the current formulation of thermodynamics for black holes is: what is the heath 
bath? When working with one temperature, Th, this is by definition the bath temperature, and normally also the 
temperature of the object. This can only apply to a black hole in equilibrium with its own Hawking radiation, which 
is an unstable and thus unphysical situation; such an approach can also not deal with black holes of different size. 
Physically there is one, and only one choice for the bath: for a black hole far away from matter, the heat bath is the 
cosmic microwave background radiation, that presently has temperature Tcmb = 2.73 K (we neglect here the small 
CMB-luctuations). So the actual problem deals with a system of which the dynamics "lives" at a second temperature, 
namely Th- As for glasses, this is a two-temperature situation, in agreement with the time scale argument: Black 
holes heavier than 10~^^Mq = 10^^ g need more time to evaporate than the present age of the universe. For them the 
evaporation process, as seen by far-away observers, is so slow, that equilibration of the cosmic background radiation 
is a fast process. 

The slow evaporation processes occur at the Hawking temperature and have entropy Sbh- The entropy of the 
background radiation outside the back hole but inside the Gedanken sphere is very small. Moreover, the emitted 
radiation will immediately leave this sphere and there is no 'back reaction', so dS^^ = 0. This simplifies eq. (||) to 

dQ = TndSBH (23) 

The standard first law dU = dQ +dW thus indeed reproduces (|2l|) , notwithstanding its non-equilibrium thermody- 
namic nature with Tbh Tcmb- According to the second law (P2[), Sbh has to satisfy 

{Tcmb - Th) dSnH > (24) 

Hawking radiation leads to dSBH < 0. Eq. (|2j) is thus satisfied as long as Th > Tcmb, but not below that. One 
might think that Tcmb plays no physical role whatsoever, and only shows up as determinator in the second law. 
However, the real point is that we not yet considered absorption of background radiation by the black hole. The 
absorption rate will be proportional to the area times the energy density, i.e., ^ M'^T^^j^- The quantum absorption 
process is simply the time-reversed evaporation process. For non-rotating, neutral holes Hawking radiation leads to a 
mass loss 

The dimensionless constant depends on the type of particles present, and their absorption probabilities, called 
"oscillator strengths" in solid state systems. Th enters through the Bose-Einstein occupation numbers (for bosons, 
in particular photons) or Fermi-Dirac occupation numbers (for fermions). For an uncharged, non-rotating black hole 
Page finds a = 5.246 x 10~^ in the high-frequency limit, and 0.181 x 10~^ in the low frequency limit For 
absorption by the black hole of a photon (or a particle) from the cosmic background, the time-reversed problem 
shows up. It thus holds that aabsiT) — aem(T), no matter the character of the particle content; for simplicity we 
shall now replace both by a constant. The only difference between the two situations is the temperature occurring 
in the occupation numbers: for Hawking emission it is Th, while for cosmic background absorption it is Tcmb- The 
combined processes of Hawking emission and background photon absorption thus yields for a neutral, non-rotation 
black hole H 

^^h-^cmb) (^0) 

It exhibits an instability at Th = Tcmb, related to the fact that the "Bekenstein-Hawking" specific heat Cbh = 
dU/dTn is negative. If there is equilibrium, and Tcmb is changed a little, then Th branches away from it. 

There are two regimes. In the "classical" regime Th < Tcmb the black hole absorbs more radiation than is emits. 
Its entropy will increase, anddQ = TndSBH > 0, but this is still in accord with the second law (p4|). In the "quantum" 
regime Th > Tcmb the black hole emits more than it absorbs. Now it holds that dSnH < 0, confirming again that 
heat flows from high to low temperature. 

In analogy with glasses, one can define the apparent specific heat C = dU/dTcMB = U /Tcmb- For black holes 
this object is less natural because the background temperature cannot be changed by hand. However, C does have a 
meaning in our expanding universe. Due to the decrease of Tcmb, there will be less and less background energy to 
be absorbed. Eq. ( p6| ) tells us that, provided Tcmb decays quicker than t~^/^, a black hole will reach its maximal 
size at some moment t — to where the temperatures match, Th = Tcmb = To; from then on it will shrink, and its 
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Hawking temperature rises. In our matter dominated Universe one has Tcmb t~^^^, so finally the black hole will 
evaporate. Around to the apparent specific heat takes a form independent of Tcmb, viz. C = ksit — to)/T, with 
characteristic time scale the quantum time r = ?i/[(167r)^aA:B2o]. In the classical regime {t < to) C is negative, while 
in the quantum regime it is positive. 

The third law of black hole dynamics is related to extremal black holes, that have Th ~ 0. The third law of 
thermodynamics, however, concerns the vanishing of the entropy for Tcmb ~^ 0. We have seen already that finally all 
black holes evaporate, thereby lowering their configurational entropy very much, in accord with Planck's third law. 
Notice that this has nothing to do with the third law of black hole dynamics. What happens ultimately with the 
black hole has been the focus of studies by 't Hooft . 

The entropy change of the universe is found as for black body radiation 

dSm _ TndSm _ TndSm _ '^TniTfj — Tqmb) j-2y-j 



dSBH dU dUm HTfj-T^j^js) 

yielding the entropy production S = Sm + Sbh 

A_ akl {Tjj + 2ThTcmb + iT^MB)iTH - Tcmb) 



2inh Til 



(28) 



We now consider the whole universe as our system, so the entropy of the universe 5*^ has to be taken into account. 
While S — Srn + Sbh is the total entropy, eq. (^) becomes dQ = TcAiBdSm + TndSBH- As dQ = 0, the second law 
again leads to (p4|), but with different entropy production. For the present case we find the new result 

g ^ «fc| (r^ " Tcmb){Th - Tcmb) 

8TTh T%TcMB ^ ' 

This expression exceeds eq. (p8|), which referred to radiation that was still located near the balck hole. The difference 
is due to loss of information on the emitted radiation in our negatively curved Universe. 

Our approach thus incorporates the known properties of dynamics, and shows how the generalized second law 
comes into the play. Both the formation and evaporation of black holes leads to an increase of the entropy of the 
whole universe. Our picture involves the standard zero-entropy formulation of the third law of thermodynamics, thus 
putting aside the third law of black hole mechanics. To the best of our knowledge, there is no contradiction with the 
occurrence of negative specific heats. 



V. GLOBULAR STAR CLUSTERS 



The success of the applying thermodynamics to black holes has made us optimistic about the possibility to do the 
same for other gravitating systems. Let us therefore discuss the case of globular star clusters. In our Galaxy some 
150 of them have been observed, while their total number is estimated to be some hundreds. Each has typically 
N = 10^-10^ stars of solar mass. For a star cluster the evaporation time Tcvao ~ lOOrioiax, with Trdax ~ -^"'^/^rcross 
is much larger than the time Tcross one star needs to cross the cluster |30) [0. This implies that thermodynamics 
involves at least two temperatures, the cosmic background temperature, and the (space dependent) local temperature 
related to the average kinetic energy of the stars. Connected to each of these temperatures there is an entropy. 

Attempts to formulate thermodynamics for star clusters have been partly successful. The start was made by Antonov 
in the early 1960's, and extended by Lynden-Bell in the late 1960's |33j. They considered a cluster located in a 
rigid sphere, and kept at fixed "bath" temperature, the so-called isothermal sphere, and studied the partition sum of 
such a system. For a recent review, see Q . Though this setup does not answer the principle problems of the instability 
of gravitating A'^-body systems, the gravothermal catastrophe was used to explain the energy distribution function 
(Lynden-BcU distribution), and describes the non-stationary ("violent") relaxation. Later it was shown that small 
perturbations of the phase trajectories of the stellar system deviate exponentially in time, which leads to a mixing 
character ("AT- mixing"). This is expressed by the positivity of the Kolmogorov- Sinai entropy of spherical A^-body 
gravitating systems and implies exponential damping of fluctuations and hence exponential rate of relaxation. |3C| ] 

For applying our knowledge on two-temperature thermodynamics to star clusters, let us consider a spherical cluster, 
with typical radius a. We first have to say what is the system to be considered. We consider a "Gedanken-sphere" 
with radius i?, with its center overlapping with the one of the cluster. The inner part of the sphere is our "system" , 
and we are only interested in what happens with this part. The sphere is not a physical object, and has no influence on 
the motion of the stars. If no work is done, then the only change can be a change of heat, dQ, through the boundary. 
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For trying to apply thermodynamics, we notice that the first law (|l|) is, of course, satisfied, as is used in all theoretical 
approaches. We thus only have to interpret the second law (||). The left hand side, dQ, has already been identified. 
Next question is: what is T ? Our Gedanken sphere has the cosmic background temperature Tcmb ^ iK. Since 
the kinetic energy of a star of one solar mass moving within a star cluster at a speed of, say, 10 km/s is equivalent 
to a "temperature" of some lO^^if, we may set the bath temperature T equal to zero. Physically this just says that 
cosmic background radiation has no influence on the motion of the stars. 

It is known that star clusters, being in quasi-equilibrium, have to evaporate. During one relaxation time a fraction 
of about one percent of the stars leaves the cluster, never to return. It is the main purpose of this section to point out 
that this effect can be viewed as an immediate consequence of the second law. The derivation is hardly a derivation, 
it is just a bold, successful application of the second law to this situation. Indeed, the second law says that heat must 
flow from high temperatures to lower ones, so in the outward direction. An other way to see this is to insert T = 
eq. (||), implying again that dQ < 0. At distances far from the cluster, R a, re-shuffling of energy between the 
stars cannot occur, so a dQ < must indeed describe evaporation. 

When one takes as system the whole universe with the star cluster in it, the system is closed, so it holds that 
dQ = 0. The second law then requires that the total entropy increases. This indeed happens, as the entropy of the 
iVosc escaped stars, essentially given by A'osc ln(_R'^/a'^), is much larger than the entropy loss of the cluster. 

We can also take the Gedanken-sphere in the interior of the cluster. For < R < a the second law still implies 
dQ < 0, with a strict inequality sign because there is no equilibrium. This describes that heat must flow from the 
inside r < R towards the outside. It can consist of two parts: first, of stars that move, on the average outwards, and, 
second, of increase in kinetic or potential energy of the outer region via energy transfer at stellar encounters. This 
heat is generated by the gravothermal catastrophe: the energy of the central region goes down, while its temperature 
goes up due to the negative heat capacity. The second law indeed requires that property since the outward energy 
flow must go hand in hand with flow from high temperatures to lower ones. Thus also the gravo-thermal catastrophy 
can be seen as an immediate consequence of the second law. 

Let us next consider a globular star cluster with an energy production in its center. A physical realization of this 
occurs when a binary star system is formed there. Indeed, a binary system can release much energy by becoming 
stronger bound. In this case we can take our Gedanken-sphere around the central region including the binary. In 
agreement with the second law, the central part looses energy, which goes to the non-central part of the cluster. As 
a result, it will become higher and higher in energy. Even if the star cluster was bound initially (i.e. having negative 
total energy) there is no reason why it must remain bound: all binding energy can be balanced by energy released 
from the binary. In short, after the accidental formation of a central binary, the whole cluster will evaporate quicker. 
This happens relatively fast, as the time scale for energy transfer to the binary is small compared to the relaxation 
time. 

When we look at shorter timescales where evaporation does not occur, we can assume thermodynamic equilibrium. 
We can then indeed insert dQ = r2dS'2 in the first law, and re-derive the standard results for isothermal spheres. For 
realistic globular star clusters one expects that a statistical description will only work if too far collapsed states are 
suppresed. In the statistical mechanics description indeed a "centrophobic force" enters, which is a van der Waals-type 
of force, that repels stars from the center. It expresses that at a finite moment in time states with too much collaps 
are unattainable, and should be suppressed in a statistical description. Gurzadyan and Keckek have considered 
thermodynamics of isolated globulars before. Starting at the level of the chemical potential, rather than from a 
partition sum, they mapped the pre-coUaps situation onto a Thomas-Fermi equation. We may arrive at the same, or 
at a polytropic equation of state. Further work on this issue is in progress. 

In conclusion, we have pointed out that the second law of thermodynamics is capable to explain the evaporation and 
the gravo-thermal catastrophy of globular star clusters, just as it explains the thermal behavior of black holes, but also 
of evaporating atmospheres and cups of water. This probably implies that thermodynamics works for self-gravitating 
systems in general, despite of the long range forces, the non-extensive behavior, and the negative heat capacities. As 
is commonly assumed, on relatively short timescales the equilibrium formulation holds. On long timescales, however, 
the non-equilibrium formulation of thermodynamics must be applied. 
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